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Abstract. An explicit derivation of dispersion relations and spectra for peri- 
odic Schrodinger operators on carbon nano-structures (including graphene and 
all types of single- wall nano-tubes) is provided. 



1. Introduction 

Carbon nano-structures, in particular fullerenes (buckyballs) , carbon nano- 
tubes, and graphene have attracted a lot of attention recently, due to their 
peculiar properties and existing or expected applications (e.g., [2D [25], 155]). 
Such structures have in particular been modelled by quantum networks (e.g., 
[21 El [261 12Z]), also called quantum graphs, which goes back to quantum graph 
models in chemistry [511 EI] and physics [U El HH HH H5] (see also [H [321 [33] and 
references therein). A direct and inverse spectral study of Schrodinger operators 
on zig-zag carbon nano-tubes was conducted in [261 EZj- 

In this paper, we take a different from [201 [2Z] approach to such a study. 
Namely, we provide a simple explicit derivation of the dispersion relations for 
Schrodinger operators on the graphene structure, which in turn implies the struc- 
ture of the spectrum and density of states. This derivation was triggered by the 
one done in [36] for the photonic crystal well by [21 [3] , albeit the pre- 

sented computation is simpler and more convenient for our purpose than the one 
in [36]. It reflects the known idea (e.g., [H [5j EU [351 EH]) that spectral analysis of 
quantum graph Hamiltonians (at least on graphs with all edges of equal lengths) 
splits into two essentially unrelated parts: analysis on a single edge, and then 
spectral analysis on the combinatorial graph, the former being independent on 
the graph structure, and the latter independent on the potential. 

The results are formulated in terms of the monodromy matrix (or rather its 
trace, also called the Hill discriminant [17j ) of the ID potential on one edge of 
the graphene lattice. Then this dispersion relation, just by simple restriction 
procedure, gives answers for any carbon nano-tube: zig-zag, armchair, or chiral. 
We would like to emphasize that relations of properties of the discriminant to the 
properties of the one-dimensional periodic potential have been studied for a long 
time and are understood by now extremely well (e.g., [121 \U\ [TBI [HI [2H El 231 
1HE31 EE] an d references therein). Thus, one can extract all spectral information 
that might be needed, from our explicit description of dispersion relations that 
involve the discriminant. 
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The paper is structured as follows: in the next Section [21 the main geometries 
and operators are introduced. Section [3] is devoted to derivation of the dispersion 
relation and the spectral structure for graphene, the main results provided in 
Theorem 13.61 The following Section H] deals with nano-tubes. The main results 
accumulated in Theorem 14.31 provide dispersion relations and all parts of the 
spectra of the nano-tube operators. The last sections contain additional remarks 
and results, and acknowledgments. 

2. SCHRODINGER OPERATORS ON CARBON NANO-STRUCTURES 

All structures studied in this text can be introduced through the honeycomb 
graphene structure [2H [2S1 ES] • So, we start with discussing the latter. 

2.1. Graphen. It is assumed that in graphene, the carbon atoms are situated 
at the vertices of a hexagonal 2D structure G shown in Fig. [T] below. We will 




Figure 1 . The hexagonal lattice G and a fundamental domain W 
together with its set of vertices V(W) = {a, b} and set of edges 
E(W) = {f,g,h}. 

assume that all edges of G have length 1. It will be crucial for us to consider the 
following action of the group Z 2 on G: it acts by the shifts by vectors piei + P2&2, 
where (pi,P2) & %> 2 and vectors ei = (3/2, y/3/2), e 2 = (0, \/3) are shown in 
Fig. [H We choose as a fundamental domain (Wigner-Seitz cell) of this action the 
parallelogram region W shown in the picture. Here two vertices V(W) = {a, b} 
are assumed to belong to W, while the vertices b', b" and b'" belong to other shifted 
copies of the fundamental domain. Three edges f,g,h belong to W. Although 
the graph G does not have to be directed, it will be convenient for us to assign 
directions to the edges in W as shown in the picture. 

We will now equip the graph G with the structure of a quantum graph (quan- 
tum network) [71 [281 1211 [331 EI] • This requires introduction of a metric structure 
and of a differential Hamiltonian. We assume that G is naturally embedded into 
the Euclidean plane, which induces the arc length metric on G, as well as the 
identification of each edge e in G with the segment [0, 1]. Under this identifica- 
tion, the end points of an edge correspond to the points and 1 (this identification 
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is unique up to a symmetry with respect to the center of the edge, i.e., up to 
the direction of the edge). We also obtain a measure (that we will call dx), and 
the ability to differentiate functions along edges and to integrate functions on G. 
In particular, the Hilbert space L 2 (G) := ® eg £(G) L 2 (e) consisting of all square 
integrable functions on G can be naturally defined. Here E(G) denotes the set 
of edges in G. 

We now describe the graphene Hamiltonian H in L 2 (G) that will be studied 
in our paper (it has also been considered for some special potentials in [2]). 
Such operators are used for approximating the band structure of carbon nano- 
structures, as well of other compounds (e.g., (2J [31 G3], and references therein). 
Let qo(x) be an L 2 -function on the segment [0, 1]. We will assume that it is even: 

q (x) = q (l-x). (2.1) 

The evenness assumption is made not just for mathematical convenience. As the 
proposition below shows, this condition is required if one considers operators in- 
variant with respect to all symmetries of the graphene lattice. One could consider 
a directed honeycomb graph, and thus avoid the evenness condition (hence losing 
invariance of the operator), but the authors did not see any compelling physical 
reason for doing so. 

Using the fixed identification of the edges e G E(G) with [0, 1], we can pullback 
the function qo(x) to a function (potential) q(x) on G. Notice, that due to the 
evenness condition imposed on qo(x), the definition of the potential q does not 
depend on the orientations chosen along the edges. It is also easy to see that the 
following claim holds: 

Proposition 2.1. The potential q defined as above, is invariant with respect to 
the full symmetry group of the honeycomb lattice G. Moreover, all invariant 
potentials from L 2loc (G) are obtainable by this procedure. 

We skip the immediate proof of this statement. 

We can now define our Hamiltonian H. It acts along each edge e as 

Hu(x) = -^^- + q(x)u(x), (2.2) 
dx z 

where we use the shorthand notation x for the coordinate x e along the edge e. 

The domain domH of the operator H consists of the functions u that belong 
to the Sobolev space H 2 (e) on each edge e in G and satisfy the inequality 

IMlW) < 00 • ( 2 - 3 ) 

egE(G) 

They also must satisfy the so-called Neumann vertex conditions (also somewhat 
misleadingly called Kirchhoff vertex conditions) at vertices. These conditions 
require continuity of the functions at each vertex v (and thus on all graph G) 
and vanishing of the total flux, i.e., 

u ei (v) = u e2 (v) if ei,e 2 6 E V (G), (2.4a) 
E <(v) = ° ( 2 - 4b ) 

e£E v (G) 
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at each vertex v G V(G). Here E V {G) := { e G E(G) \ v G e } is the set of edges 
adjacent to the vertex v, u e is the restriction of a function u to the edge e, and 
u' e (v) denotes the derivative of u e along e in the direction away from the vertex 
v {outgoing direction). Thus defined operator H is unbounded and self-adjoint 
in the Hilbert space L 2 (G) [281 El]- Due to the condition on the potential and 
Proposition I2.1[ the Hamiltonian H is invariant with respect to all symmetries 
of the hexagonal lattice G, in particular with respect to the Z 2 -shifts, which will 
play a crucial role in our considerations. 

2.2. Nano-tubes. We provide here a very brief introduction to carbon nano- 
tubes. One can find more detailed discussion and classification of nano-tubes in 
ES]. We also emphasize that only single- wall nano-tubes are considered. 
Let p G 1R 2 \{0} be a vector that belongs to the lattice of translation symmetries 
of the honeycomb structure G. In other words, G + p = G. Let us denote by 
t p the equivalence relation that identifies vectors z\,z% G G such that z^ — z\ is 
an integer multiple of the vector p. A nano-tube T p is the graph obtained as the 
quotient of G with respect to this equivalence relation: 

T p := G/ip. (2.5) 

This graph is naturally isometrically embedded into the cylinder ~R 2 /l p . 

If p = p\t\ +P2&2, we will abuse notations denoting T p by T(p ltP2 y For example, 
T(o,n) is the so-called zig-zag nano-tube, while T(n,n) is the so-called armchair 
nano-tube. The names come from the shape of the boundary of a fundamental 
domain (cf. Fig. [2]). There are many other types of nano-tubes, besides the zig-zag 




Figure 2. A zig-zag (left) and armchair (right) nano-tube TJ^) 
and Tnu, respectively. The vectors show the translation vector p. 
The name-giving fundamental domain of each of the nano-tubes is 
shaded in dark grey. The dashed lines have to be identified. Below, 
the (degenerate) zig-zag nano-tube TJ^i) is shown. 



and armchair ones. They are usually called chiral. 

A degenerate example is given by the zig-zag nano-tube 7(0,1), which consists 
of a sequence of loops ("beads") joined by edges into a ID-periodic necklace 
structure (see Fig. El). 
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One can notice that due to existence of rotational and mirror symmetries of 
the hexagonal structure G, different vectors p can produce the same nano-tubes 
T p . For instance, T (m>n ) = T( n , m ). Also, zig-zag tubes T (0jiV) , T (A r i0 ) , and T( N _ N ) 
are the same (see [21] and references therein for a more detailed classification of 
nano-tubes). 

The Hamiltonian H p onT = T p is defined exactly as for the graphene lattice G. 
Alternatively, one can think of H acting on functions on G that are periodic with 
the period vector p (this definition requires some precision, since such functions 
do not belong to L 2 (G)). 

3. Spectra of graphene operators 

Here we calculate the dispersion relation and thus all parts of the spectrum 
of the graphene Hamiltonian H (see also [2J [3] for related considerations). One 
can notice that the density of states is determined by the dispersion relation, and 
thus when the latter is known, the former can be determined as well (U 153"] . 

We apply now the standard Floquet-Bloch theory (e.g., pZ3 EU [53]) with re- 
spect to the Z 2 -action that we specified before. This theory also holds in the 
quantum graph case, (see, e.g., [El [151 [301 [351 09] an d references therein). This 
reduces the study of the Hamiltonian H to the study of the family of Bloch 
Hamiltonians H 9 acting in L 2 (W) for the values of the quasimomentum 9 in the 
Brillouin zone [— 7r,7r] 2 . Here the Bloch Hamiltonian H 9 acts the same way H 
does, but it is applied to a different space of functions. Each function u = {u e } in 
the domain of H 9 must belong to the Sobolev space u e G H 2 (e) on each edge e and 
satisfy the vertex conditions (j2.4p . as well as the cyclic conditions (Floquet-Bloch 
conditions) 

u{x + piei + p 2 e 2 ) = e ip - 9 u{x) = e iipi9l+P2(>2) u(x) (3.1) 

for any vector p = (pijPz) G Z 2 and any x G G. 

Due to the conditions (13. ip . functions u are uniquely determined by their re- 
strictions to the fundamental domain W . Then conditions (12 .4p and (13.11) reduce 
to 

'u f (0) = u g (0) = u h (0) =: A 
u' f (0) + <(0) + <(0) = 

u f {l)=e^u g {l)=e^u h {l)=:B ' K '> 
u' f (l) + e^u' g (l) + e l9l u' h (l) = 0. 

By standard arguments (see e.g. [3H Theorem 18]), H 9 has purely discrete spec- 
trum o-(H 9 ) = {Xj(0)}. The graph of the multiple valued function 9 \— > {\j(9)} is 
known as the dispersion relation, or Bloch variety of the operator H. It is known 
[T71 12~4"l I3"lj 153"! that the range of this function is the spectrum of H: 

a(H)= |J a(H 9 ). (3.3) 

6»e[-7r,7r] 2 

It is also well known that the dispersion relation determines not only the spec- 
trum, but the density of states of H as well [H [53] . 

So, our goal now is the determination of the spectrum of H 9 and thus the 
dispersion relation of H . In order to determine this spectrum, we have to solve 
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the eigenvalue problem 

H e u = Xu (3.4) 

for A G M and a non-trivial function u G L 2 (W) satisfying the above boundary 
conditions. 

Let us denote by S D the spectrum of the Dirichlet Hamiltonian if D acting as 
in 02.21) on (0, 1) with boundary conditions w(0) = u(l) = 0. If A ^ £ D , there 
exist two linearly independent solutions <po, (p\ (depending on A) of the equation 

- if" + q <p = Xip (3.5) 

on (0, 1), such that y?o(0) = 1, y3 (l) — and <^i(0) = 0, <^i(l) = 1- F° r example, 
if go = and A > then we have A ^ S D if and only if /z := a/A ^ 7rZ and 

sinu(l — t) , sinut . 

(p (t) = ^ '- and <px(t) = — . (3.6) 

sin jj, sin [i 

We will often use the notation <^ ,a, V?i,a to emphasize dependence on the spectral 
parameter. 

We will assume that the functions tpj are lifted to each of the edges in W, 
using the described before identifications of these edges with the segment [0, 1]. 
Abusing notations, we will use the same names tpj for the lifted functions. Then, 
for A ^ TP we can use (13.21) to represent any solution u of (13.41) from the domain 
of H e on each edge in W as follows: 



u g = Atp + e- iei 5^i (3.7) 
u h = Aipo + e-'^Bif!. 

With this choice, the continuity conditions in ( 13.21) and the eigenvalue equation 
on each edge are automatically satisfied. The remaining two equations that guar- 
antee zero fluxes at the vertices, lead to the system 



(3.8) 



&p' (0)A + (1 + e~ iei + e- W2 )p[(0)B = 
1 + e *i + e ie2 )if' {l)A + 3(p[(l)B = 0. 

Using the symmetry ( 12.11) of the potential qo, we obtain 

¥>i(l) = -Vo(0) and ^(l) = Vi(0). (3.9) 
Thus, the system ( 13.81) reduces to 

-3(p[(i)A + (1 + e~ idl + e- ie2 )cp[(0)B = 
[1 + e Wl + e ie2 )ip[(0)A - 3^(1)5 = 0. 

The quotient 



(3.10) 



,(A)-^-^ (3 11) 

is well defined, since <p'i(0) ^ 0. Thus, the system (13.101) can be rewritten as 

-3r](X)A + (1 + e~^ + e~^)B = 
l + e Wl +e W2 )A-3ri(X)B = 0. 
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The determinant of this system is clearly equal to 

|l + e i91 +e ie2 | 2 - V(A) 
= 3 + 2 cos 9 X + 2 cos 9 2 + 2 cos (9 1 - 9 2 ) - 9r/ 2 (A) (3.13) 
= 1 + 8 cos ^rr^ cos ^ cos f - 9?7 2 (A). 

Thus, we have proven the following lemma: 

Lemma 3.1. If X ^ S D , then X is in the spectrum of the graphene Hamiltonian 
H if and only if there is 9 G [— tt, tt] 2 such that 

9r] 2 (X) = \l + e Wl +e i92 | 2 , 

or 

a n n n 

9rj 2 (A) = 1 + 8 cos — cos — cos — . (3.14) 

Remark 3.2. Note that this lemma gives a nice relation between the spectrum 
of the metric graph Hamiltonian H and the discrete graph Laplacian [9l [10] 
AU := dk^J2 w ~v U ( w )- Indeed, the Bloch Laplacian A 9 on £ 2 ({a,b}) = C 2 has 
the matrix 

If 1 + e- Wl + 



A 6 

3 V 1 + e ' ei + eie2 
Thus, the lemma is just the statement that for A ^ S D , we have A G o~{H e ) if 
and only if 97(A) G a(A e ), and hence A G cr(H) if and only if 77(A) G cr(A). 

This relation between quantum and combinatorial graph operators is well 
known and has been exploited many times (e.g., [H El [H [3U [351 E3 15D]). 



Lemma 13. 1[ in particular, says that in order to find the spectrum of H, we 
need to calculate the range of the following function on [— tt, tt} 2 : 



F{0i, 9 2 ) = ^l + 8 cos 9j -^ 1 cos y cos y- ( 3 - 15 ) 

Lemma 3.3. T/ie function F has range [0,3]; its maximum is attained at (0,0) 
and minimum at (2tt/3, — 27r/3) and (— 27r/3, 2tt/3). 

The proof of the lemma is straightforward, after noticing that the function 
F(9 h 9 2 ) 2 = 1 + 8 cos cos | cos | = |1 + + e^\ 2 

ranges from to 9. 

Next, we want to interpret the function 77(A) in terms of the original potential 
qo(x) on [0, 1]. To this end, let us extend qo(x) periodically to the whole real axis 
ffi. and consider the Hill operator H peT on M given as in (12.21) with the resulting 
periodic potential: 

H**u(x) = + qo(x)u(x). (3.16) 

(We maintain the notation qo(x) for the extended potential.) 
We will be interested in the well studied spectral problem 

H pev u = Xu (3.17) 
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As it is usually done in the theory of periodic Hill operators (e.g., [T7J [2U SH 
021 E3]), we consider the so called monodromy matrix M(A) of H per . It is defined 
as follows: 




where <p is any solution of (13.171) . In other words, the monodromy matrix shifts by 
the period along the solutions of (13.171) . The matrix valued function A i— > M(A) 
is entire. 

It is well known (see the references above) and easy to see that the mon- 
odromy matrix has determinant equal to 1. Its trace plays the major role in 
the spectral theory of the Hill operator. Namely, let -D(A) = trM(A) be the so 
called discriminant (or Lyapunov function) of the Hill operator H pei . The next 
proposition collects well known results concerning the spectra of Hill operators 

PH EHUD E2l E3]: 

Proposition 3.4. 

(i) The spectrum a(H per ) of H per is purely absolutely continuous. 

(ii) a{H pev ) = { A e R | |D(A)| < 2 }. 

(hi) a(H per ) consists of the union of closed non-overlapping and non-zero 
lengths finite intervals (bands) B 2k ■= [a 2 k,b 2 k\, B 2k +i ■= [&2fc+i, a 2k+1 ] 
such that 

ao < bo < bi < ai < a 2 < b 2 < ■ ■ ■ 
and linn^oo a k = oo. 

The (possibly empty) segments (b 2 k,b 2 k+i) and (a 2 k+i, 02^+2) o,re called 
the spectral gaps. 

Here {a k } and {bk} are the spectra of the operator with periodic and 
anti-periodic conditions on [0, 1] correspondingly. 

(iv) Let A]? G S D be the kth Dirichlet eigenvalue, labelled in increasing order. 
Then A^ belongs to (the closure of) the kth ga$\. When go is symmetric, 
as in our case, A^ coincides with an edge of the kth gapj. 

(v) // A is inside the kth band Bk, then -D'(A) 7^ ; and D(\) is a homeo- 
morphism of the band Bk onto [—2,2]. Moreover, D(X) is decaying on 
(—00, bo) and (a 2 k, b 2 k) and is increasing on (62/c+i, «2fc+i)- It has a simple 
extremum in each spectral gap [0^,0^+1] an d [bk,bk+i]. 

(vi) The dispersion relation for H peT is given by 

D(\) = 2cos6, (3.18) 
where 9 is the one- dimensional quasimomentum. 

There are many other important direct and inverse spectral results concerning 
the well studied (in particular, due to the inverse scattering transform research) 
operator H peT (see, e.g., [T71 [111 [IHl [211 Sll SI1 S31 E31 [58] ) . 

We will now see the relation between the function f](X) that was introduced 
for the graphene operator H and the discriminant of H pcT . In order to do so, we 
introduce another basis of solutions of (13.171) . namely c\ and s\ with c\(Q) = 1, 



If the gap closes, this boils down to a single point. 
The same comment applies here. 
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c' A (0) = and s x (0) = 0, s' x (0) = 1 (the notations are chosen to remind the cosine 
and sine functions in the case of zero potential). Using this basis of the solution 
space, we obtain 




(3.19) 



A simple calculation (assuming the symmetry (12. ip ) relates the new basis with 
the one of ip and <p\ (which we now denote ^oA'V^i.a to emphasize dependence 
on the spectral parameter): 

c\ = <A),a + »7(A)y»i,A and s x = , fr .. (pi,\- 

In particular, c\(l) = s' x (l) = //(A). Thus, 

17(A) = ~L>(A). (3.20) 

For example, if go = 0, then 77(A) = cos(vT). 

So far, we have been avoiding points of the Dirichlet spectrum S D of a single 
edge (i.e, the Dirichlet spectrum of the potential qo(x) on [0,1]). We will now 
deal with exactly these points. 

Lemma 3.5. Each point A G TP is an eigenvalue of infinite multiplicity of the 
graphene Hamiltonian H . The corresponding eigenspace is generated by simple 
loop states, i.e., by eigenfunctions which live on a single hexagon and vanish at 
the vertices (see Fig. below). 

Proof. We need to guarantee first that each A G S D is an eigenvalue. Indeed, an 
eigenfunction is provided by a simple loop state of the type shown for zero poten- 
tial in Fig. [3] below. It is constructed as follows. Let ip\(x) be an eigenfunction 



A! 




Figure 3. A simple loop state constructed from an odd eigenfunc- 
tion on [0, 1]. 

of the operator —d 2 /dx 2 + qo(x) with the eigenvalue A and Dirichlet boundary 
conditions on [0, 1]. Then, due to the symmetry (evenness) of the potential, we 
can assume the eigenfunction to be either even, or odd. For an odd eigenfunc- 
tion ip x (x), repeating it on each of the six edges of a hexagon, we clearly get 
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an eigenfunction of the operator H. If it is an even eigenfunction, then repeat- 
ing it around the hexagon with an alternating sign does the same trick. Thus 
A G <Tpp(H). 

It is well known then (e.g., |17j ) that the multiplicity of the eigenvalue must be 
infinite. For completeness, we repeat here this simple argument. Let M\ C L 2 (G) 
be the eigenspace. Consider a vector 7 that is a period of the lattice G and the 
operator 5 7 of shift by 7 in L 2 (G). Then S 7 acts in M\ as a unitary operator. 
If dim Ma were finite, S 1 would have had an eigenvector / G M\ C L 2 (G) with 
an eigenvalue /1, = 1. However, such a function / obviously cannot belong 
to L 2 (G), since it is quasi-periodic in the direction of the vector 7. This proves 
infinite multiplicity. 

We note now that due to [30| (see also [33 Thm. 11]), linear combinations 
of compactly supported eigenfunctions are dense in the whole eigenspace M\. 
Thus, we only need to show that the simple loop states just described generate 
all compactly supported eigenfunctions in the space M\. 

Suppose that cp is a compactly supported eigenfunction of H corresponding to 
the eigenvalue A G TP. First, note that ip vanishes at each vertex. Indeed, if this 
were not true, due to compactness of support, there would have been an edge 
such that at its one end v (corresponding to x — 0), (f(v ) 7^ 0, while at the 
other V\ (corresponding to x — 1), <p(vi) = 0. Expanding into the basis c\, s\, we 
get 



In particular, <p(0) = A ^ 0. Then <p(l) = Ac x (l) = As' x (l) ^ 0, since s A (l) = 0. 
This leads to a contradiction. Thus, <p vanishes at all vertices. In particular, 
on each edge it constitutes a Dirichlet eigenfunction for the Hill operator on this 
edge. 

This also implies that the support of ip, as a graph, cannot have vertices of 
degree 1 and supp ip cannot be a tree. Thus, there must be a loop in the support 
of ip. In particular, the outer boundary of the support must be a loop. Take one 
boundary edge Cq. There is a hexagon inside the boundary loop which contains 
this edge. Consider a simple loop state (p coinciding with the eigenfunction ip 
on e and extended to the hexagon the way it was described before. Subtracting 
(p from ip, we obtain a new eigenfunction <p. The number of hexagons inside the 
boundary loop of the support of <p is less than inside the support of <p. Thus, 
continuing this procedure, we eventually represent <p as a combination of simple 
loop eigenstates. Fig. @] below illustrates this process. □ 



<p(x) = Ac\(x) + Bs x (x). 




Figure 4. Deleting simple loop states (dark grey) from the sup- 
port of an eigenfunction (light grey). 
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We can now completely describe the spectral structure of the graphene operator 
H. 

Theorem 3.6. 

(i) The singular continuous spectrum o~ sc (H) is empty. 

(ii) The absolutely continuous spectrum a ac (H) has band-gap structure and 
coincides as a set with the spectrum o~(H pcv ) of the Hill operator H per 
with potential q periodically extended from [0, 1]. In particular, 

a ac (H) = {XeR\\D(X)\<2}, 

where D(X) is the discriminant of H peT . 

(iii) The pure point spectrum a pp (H) coincides with TP and thus, due to the 
evenness of the potential, belongs to the union of the edges of spectral 
gapsofa(H^) = a ac (H). 

(iv) The dispersion relation consists of the variety 

„ / . n 2 / (9i — 6o 6i 8o , . 

D A = ±-\l + 8 cos- cos — cos— 3.21 

v ; 3V 2 2 2 K J 

and the collection of flat branches A G £ D . 

(v) Eigenvalues A G S D of the pure point spectrum are of infinite multi- 
plicity and the corresponding eigenspaces are generated by simple loop 
eigenstates. 

In particular, a(H) has gaps if and only if a(H peT ) has gaps. 
Proof. 

The claim (i) about absence of the singular continuous spectrum is a general 
fact about periodic "elliptic" operators. For instance, the standard proof applied 
for the case of periodic Schrodinger operators in j53j [56] or [3U Theorem 4.5.9] 
works in our situation. Alternatively, in [20] one can find this statement proven 
for a general case of analytically fibered operators, which covers our situation as 
well. 

Statement (iv) is a combination of Lemmas 13.11 and 13.51 and formula (I3.20p . 

The statement (ii) about absolute continuity of the spectrum outside the points 
of S D follows from (iv) and the standard Thomas' analytic continuation argu- 
ment [311 E3 [56]. We remind the reader that according to this argument, eigen- 
values correspond to constant branches of the dispersion relation. It is clear that 
the dispersion curves f)3.2ip have no constant branches outside S D . 

The claim (iii) is just a combination of Lemma I3~51 and of the Thomas' argument 
again, which shows that there are no eigenvalues outside T D . 

Finally, (v) is a combination of (iii) and Lemma 13.51 □ 

It is clear that the function F 2 has non-degenerate minima F = at the points 
9 = ±(27r/3, — 27r/3). Thus, the function ±F has conical singularities at these 
points. Since, according to Proposition I3.4[ the discriminant D is monotonic 
with non-degenerate derivative near each point where -D(A) = 0, one obtains the 
following conclusion: 

Corollary 3.7. The dispersion curve of the graphene operator H has conical 
singularities at all spectral values A such that D{\) = 0. 
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These singularities (sometimes described as "linear spectra") represent one of 
the most interesting features of graphene structures (cf. Fig. These singu- 




FlGURE 5. The dispersion relation cos vA = ±F(6)/3 for the free 
(i.e., go = 0) case. The cones arising from D(X) = 2cos(a/A) = 
are at the levels A = (7r(2fc + 1)) 2 . They are located inside a band of 
the corresponding Hill operator. Note that the cones at -D(A) = ±2 
(i.e., A = (irk) 2 in the free case) are located at the band edges of 
the Hill operator. 

larities resemble Dirac spectra for massless fermions and thus lead to unusual 
physical properties of graphene (e.g., [25]). We see that quantum graph models 
with arbitrary periodic potentials preserve this feature. 



4. Spectra of nano-tube operators 

We use here the notations concerning nano-tubes that were introduced in Sec- 
tion [2T2J Consider a vector p = (^1,^2) £ ^ 2 that belongs to the lattice of 
translation symmetries of the graphene structure G, i.e., it shifts the structure 
by Pi + P2&2 (see Fig. [Q. We will use, as before, the corresponding nano-tube 
T p = T(p liP2 ) and the Hamiltonian H p = if( Pl P2 ) on T p (see Section I2T21 . 

Let B = [— 7r,7r] 2 be the Brillouin zone of graphene. Then, as we discussed 
in the previous section, the Floquet-Bloch theory provides the direct integral 
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expansion 

H = I H e d9. (4.1) 



B 



In the case of the nano-tube T p , only the values of quasimomenta enter that 
satisfy the condition 

p ■ 6 = + p 2 9 2 G 2vrZ (4.2) 
since a function on T p lifts to a p-periodic function u on G, i.e., 

u(x + pxe x + p 2 e 2 ) = u(x) 
(cf. ( 13. ip ). Thus, let us consider the following subset B p C B: 

B p := {6 G [-7r,7r] 2 \p-6 G 2ttZ}. (4.3) 

Then, we have the direct integral decomposition 

e 

H p = J H e d9. (4.4) 

Bp 

In particular, the spectrum of H p is given by 

a(H p ) = |J a(i7*), (4.5) 

and the dispersion relation for H p is just the restriction to B p of the dispersion 
relation for H described in the part (iv) of Theorem 13.61 

This implies for instance that we still have TP C <j pp (H p ) and the rest of the 
spectrum is determined by the pre-image 

iT^iiiW} = D^{± 2 -F{B P )). (4.6) 

One should notice that it is conceivable that non-constant branches of the dis- 
persion curves (I3.14p might sometimes have constant restrictions to B p , thus 
providing new eigenvalues for the nano-tube Hamiltonian. This happens if the 
line ( 14. 2 p is a level set of the function F. It is easy to find such lines. 

Lemma 4.1. The only linear level sets of the function F inside B are the fol- 
lowing ones: 

(i) 0j = ±7T 
(ii) 6 2 = ±7T 

(iii) X - 2 = ±7T 

On these lines F(8i,6 2 ) = 1. 

The proof of the lemma is immediate from the expression (I3.15P for the function 
F (see also Fig. [6] for the level sets of F, which illustrates this statement). We 
will see that existence of such lines leads to additional pure point spectrum for 
some types of nano-tubes. 

In order to determine the spectra of nano-tubes, we need to know the ranges of 
the function F (cf. ( 13.151) ) restricted to B p . These are described in the following 

Lemma 4.2. 
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FIGURE 6. The level curves of F for levels varying from (the two 
dots at ±(27r/3, — 27r/3)) till 3 (the dot at the origin). The level 
curve associated to F — 1 is dotted, the areas where F < 1 are 
shaded. The dashed line is the line in B p closest to the minimum 
point (— 27r/3, 27r/3) (in the case pi — p 2 = 3m ± 1). 



(i) The function F restricted to B p , achieves its maximum 3 at (0, 0) G B p 
for any p. 

(ii) The location and value of the minimum depends on the vector p. The 
minimal value 

a{p) := min F{9) 

satisfies 

«(p)e[o,i] (4.7) 

for any p. 

(hi) a(p) = if and only if pi — p 2 is divisible by 3. 

(iv) a(p) = 1 if and only if p = (0,±1), (±1,0), (0,±2), (±2,0), (1,-1), 
(—1, 1), (2, —2), or (—2,2). (All these cases correspond to zig-zag nano- 
tubes). 

(v) In all cases not covered by (Hi) and (iv), let pi — p 2 = 3m ± 1. Then 
a(p) G (0, 1) can be found by minimizing the function F over the line 

+P2&2 = 2-rrm. 

In particular, in the case when p = (0, N) with N = 3m ± 1 > 2 
(m E Z), one has 

77-/77 

a((0,AO) = 2 cos— -1 (4.8) 
(this formula gives the correct answer a(p) = also when N = 3m). 
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(vi) 



lim a(p) = 0. (4.9) 

\p\— *oo 



Proof. 

(i) The claim about the maximum is straightforward, since the only maximum 
point (0, 0) of F in B belongs to B p for any p. 

(ii) The expression (13 . 1 5[) shows that the set of points 9 G B where F — 1, 
consists of four lines 9j = ±7r, as well as two lines 9\ — 9 2 = ±7r. Since no 
line p ■ 9 = can miss all these points, we conclude that a(p) < 1 for any p. 
The inequality a(p) > is obvious, since, as we have already discovered, the 
expression under the square root in (13 . 15)) has its minimum equal to 0. 

(iii) As we have already indicated before, the points where F reaches its min- 
imum are (2tt/3, — 2n/3) and (— 27r/3, 2n/3). Thus, for a(p) = to hold, one of 
the lines p-9 = 2nn must pass through one of these points. Thus, p\—p 2 = ±3n, 
and the claim is proven. 

(iv) In order for a(p) to be equal to 1, the lines p ■ 9 G 27rZ should not enter 
the zones where F < 1 (shaded in Fig. [6]). It is clear that when pi,P2 are of the 
same sign, this is impossible for the line p ■ 9 = 0, unless one of the coordinates 
Pj is equal to zero. One can assume then, due to symmetries, that p = (0,N), 
N > 0. In this case, the line p ■ 9 = N9 2 = 2tt enters the shaded region, unless 
N < 2. 

If the coordinates pj have opposite signs, then in order for the first "non-trivial" 
line p ■ 9 = ±2ir not enter the Brillouin zone (and thus in particular the shaded 
area), one has to satisfy the condition \pin — p 2 n\ < 2n. Due to the signs of pjS 
being opposite, this means that \px\ + \p 2 \ < 2. This restricts the situation to 
the vectors p = (1, —1) and (—1, 1) only, which do satisfy a(p) = 1. If this line 
does enter the Brillouin zone, the only case when the shaded area is not entered 
is when p = (2, —2) or (—2, 2) and thus the line goes along the boundary of the 
shaded region. 

(v) In order to find ac(p), we need to minimize F over the set of lines p- 9 = 2-nn 
for all such integers n that the line intersects the Brillouin zone B. This entails 
first determining the appropriate of value of n and then minimizing over the 
corresponding line. The minima of F are located at the points (— 27r/3, 27r/3) 
and (27r/3, — 27r/3) and are the only local minima in the shaded regions shown 
in Fig. Thus, we need to find the value of n that provides a line closest to a 
minimum point (see again Fig. [H]). Evaluating p ■ 9 at the point (— 27r/3, 2tt/3), 
we get 2n{pi — p 2 )/3 = 2n(m ± 1/3). This suggests that the line p ■ 9 = 2-nm is 
the right one. 

When p — (0, N) with N = 3m ±1 > 2, we conclude that we need to minimize 
F over the line N9 2 = 2Ttm. Substituting the value 9 2 = 2Trm/N into the modified 
expression for F, 



F(0 U 9 2 ) = y 1 + 4 cos -j [(cos -j + cos(0! - £)] , (4.10) 

leads to a simple minimization with respect to 9\ and thus to the formula (14.81) . 

(vi) This claim is clear, since when \p\ —>■ oo, the lines p-9 G 27rZ form a dense 
set in B, and thus the minimum of F over B„ approaches zero. □ 
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Let us now concentrate for a moment on the additional pure point spectrum 
that arises due to the linear level sets described in Lemma 14.11 Let us assume 
that p = (0, 2N), AfeZ. Then, according to that lemma, the line p ■ 9 = 2iW is 
a level 1 set of F{6). Consider A such that 77(A) = 1/3 (or 77(A) = -1/3). We will 
now construct a compactly supported eigenfunction for H^n)- In order to do 
so, let us notice that 77(A) = 1/3 means that <p[ A (0) = 3^ A (l). Let us construct 
a function ip + on the boldface structure Z in Fig. [7J below. It is constructed as 




Figure 7. The extra eigenstates outside the Dirichlet spectrum 
for the zig-zag nano-tubes with even number of "zig-zags". On 
the left, the support is shown; on the right, the eigenfunctions tp± 
corresponding to the smallest solution of 77(A) = cosvA = ±1/3 
are plotted (in the case of zero potential qo = 0). 

follows: on the two "horns" directed toward the vertex a, we define the function 
to be equal to ipi t \(x). It is similarly defined on the "horns" leading towards b. 
On the "bridge" between a and b, we define the function as y?o,A + V^i.a- ^ i s eas y 
to conclude that the equality tp' x A (0) = 3^ A (l) implies that the function satisfies 
Neumann conditions at both vertices (and certainly the equation Fop = \ip on 
the edges of Z). The graph of the function <p + is visualized in the middle of 

Fig. m 

Analogously, if 77(A) = —1/3, one creates a function ip_, changing the value on 
the bridge to y?o,A — <fi,x (see the right graph in Fig. [7]). 

The functions (p± are defined on Z only, but can be extended to the whole 
structure G as follows: one repeats the functions up and down (to dashed-line 
hexagons in Fig. [7J, alternating the sign. Outside this column of hexagons, we 
define the functions to be equal to zero. These functions are periodic with respect 
to the vector e 2 with period 2, and thus define compactly supported eigenfunc- 
tions for any even zig-zag nano-tube Tr 0j2N y We will call such eigenfunctions the 
three-leaf eigenfunctions (the name suggested by the right graph in Fig. [7J). 

We are now ready to establish the main result about the spectra of carbon 
nano-tubes. First of all, let us collect all the notions we need here. As before, 
P = (Pi)P2) € Z 2 is a translation vector that determines the nano-tube T p . The 
Hamiltonian H p on L 2 (T P ) is defined as before, using the pull-back of a potential 
qo(x) on [0, 1] symmetric with respect to the point 1/2. We also denote by D(\) 
the Hill discriminant (trace of the monodromy matrix) of the Hill operator H pcr 
on R with periodized potential q Q . The subset B p of the Brillouin zone B is 
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defined in ( 14.31) . Finally, the function F(6) of the quasimomentum 9 is defined 
in (I3.15p . and a(p) is described in Lemma [4.21 

In order to avoid lengthy formulation, in the Theorem and two Corollaries 
below, when dispersion relations are described, the flat branches corresponding 
to the pure point spectrum are omitted. The pure point spectra (and thus the 
flat branches) are described in separate statements. 



Theorem 4.3. 

(i) The non-constant part of the dispersion relation for H p is provided by 

2 
3 

(ii) The singular continuous spectrum a sc (H p ) is empty. 

(iii) The absolutely continuous spectrum is given by 

2 , xn r 2 



D(X) = ±-F(6), 6 e B p . (4.11) 



aUH p ) = D- 1 ([-2, --«(?)] U [-a(p),2]) (4.12) 



and 



^([-2,-1 U [|,2]) C a ac (H P ) C a(H^) = D' 1 ([-2, 2]). (4.13) 

(iv) cr ac (ff p ) = a(H pcr ) if and only if pi — pi is divisible by 3. 

(v) a ac (H p ) = .D _1 ([— 2, — |] U [|,2]) if and only if T p is either a (0, l)- ; or 
a (0,2)- zig-zag nano-tube (or equivalent, e.g. Tn^y). 

(vi) Unless T p is a zig-zag nano-tube with an even number of zig-zags (i.e., 
T( ,2N)), one has 

a pp (H p ) = S D . (4.14) 

This spectrum consists of one edge of each spectral gap (including the 
closed ones) ofa(H per ). All these eigenvalues are of infinite multiplic- 
ity and the corresponding eigenspaces are spanned by simple loop eigen- 
functions (supported on a single hexagon) and tube loop eigenf unctions 
(supported on a loop around the tube). 

(vii) IfT p is a zig-zag nano-tube with an even number of zig-zags (i.e., T^^n)), 
one has 

a pp (H p )=Z D U~, (4.15) 

where 

E = D^ 1 ( 4 

3, 

The eigenvalues from E are of infinite multiplicity, are embedded into 
o- ac (H p ), and are located two per each band ofa(H per ). The corresponding 
eigenspaces are generated by the compactly supported three- leaf functions. 

The eigenvalues from TP are of infinite multiplicity and the corre- 
sponding eigenspaces are spanned by simple loop eigenfunctions and tube 
loop eigenfunctions. 

(viii) lfp\—pi is divisible by 3 ; the ac-spectrum of H p has exactly the same gaps 
as a(H per ). Otherwise, there is an additional gap D _1 ((— |a(p), |a(p))) 
inside each band of the spectrum of H per . 



D- 1 (±-). (4.16) 
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Proof. 

The first statement coincides with f)4.6p . 

The claim (ii) is proven exactly as the corresponding statement in Theorem l3.6l 
Statements (iii) through (v) follow from (i) and Lemma 14.21 
The statement (vi) is almost completely proven, except the description of the 
eigenfunctions. The proof of this description works exactly like in Theorem 13. 6^ 
except that the procedure of eliminating hexagons does not have to end with an 
empty set. One can end up with a loop of edges around the tube, which thus 
does not encircle any hexagons. This would provide a loop eigenfunction that 
runs around the tube, rather than around a hexagon. The similar claim in (vii) 
concerning the eigenvalues from E D is proven exactly the same way. 

What remains to be proven in (vii), is the structure of the eigenfunctions cor- 
responding to A G H. It is proven similarly to the structure of eigenfunctions 
corresponding to TP. Indeed, again according to [35J, the eigenspace is spanned 
by compactly supported eigenfunctions. Consider the outer boundary of the sup- 
port and start eliminating hexagons inside as follows. There must be a vertex 
(like the ends of horns in Fig. [7]) that borders zero values. Then, on the corre- 
sponding horn the function must be proportional to (pi : \. Let us now extend it 
to a three-leaf eigenfunction and subtract from the original one. Continuing this 
process, we eventually eliminate all hexagons. Notice that in this case we cannot 
end up with a loop around the tube, since this would force all the vertex values 
to be equal to zero, which is impossible, when A does not belong to S D . Thus, 
only the three-leaf states enter the eigenfunction. 

The statement (viii) follows from the previous ones. □ 

We will specify this result for the cases of zig-zag (p = (0, N)) and armchair 
(p = (N, N)) nano-tubes. The zig-zag case was also considered in 



Corollary 4.4. Let T^at) be a zig-zag nano-tube with N zig-zags. 

(i) The non-constant part of the dispersion relation for H(Qm is given by 



D(X) = ±3^1 + 4cos — [(cos — + cos (5i - — J J , (4.17) 

where < n < N. 

(ii) The singular continuous spectrum is empty. 

(iii) The absolutely continuous spectrum is given by 

<r aB (H (0tN) )=D- 1 ([-2,~a] U [~a,2]) (4.18) 

where a = a((0,N)) e [0,1] is defined in (14.81) . In particular, a = (i.e., 
o~a,c(H(o,N)) = c(H pcr )) if and only if N is divisible by 3. Furthermore, 
a = 1 if and only if N = 1 or N = 2 . 

(iv) If N is odd, then the pure point spectrum is given by a pp (H^ ^) = E D . 
The corresponding eigenspaces are infinite-dimensional and generated by 
simple loop eigenfunctions (supported on a single hexagon) and tube loop 
eigenfunctions (supported on a loop around the tube). 

If N is even, then a pp (H^^)) = S D U S where H is defined in (14.161) . 
In particular, if N = 2 then the embedded eigenvalues from 5 are located 
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at the band edges of a^Hp). If N > 2 is even, this new point spec- 
trum is located inside the bands. The eigenspaces corresponding to 5 are 
infinite- dimensional and generated by the compactly supported three-leaf 
functions. 

(v) a(-?f(o,7v)) has additional gaps (other than the gaps of a{H pcr ) ) if and only 
if N is not divisible by 3. 

Remark 4.5. In order to avoid confusion, we need to specify what a simple loop 
eigenstate is for the case of the necklace tube 7)o,i)- I n this case, the image of 
a hexagon in the tube is a "dumbbell" consisting of two beads of the necklace 
connected with a link (see Fig. [8] below). A simple loop eigenfunction in this case 



Figure 8. A dumbbell image of a hexagon. 

can concentrate either on a single bead, or on the whole dumbbell. 

Corollary 4.6. Let T^ N>N ) be an armchair nano-tube. 

(i) The non- constant part of the dispersion relation for H^n^n) is given by 

D(X) = ± ^l + 8cos( 9l -^)co S (|)co S (|-^), (4.19) 

where < n < N . 

(ii) The singular continuous spectrum is empty. 

(iii) The absolutely continuous spectrum is given by 

^ c (H(n,n)) = J D- 1 ([-2,2]) = a(H^). (4.20) 

(iv) The pure point spectrum is given by a pp (H^^^) = TP and is located at 
an edge of each gap in c^H^ ^) . The eigenvalues are of infinite mul- 
tiplicity and the eigenspaces are generated by simple loop eigenf unctions 
(either on a single hexagon or a loop around the tube). 

(v) a(H(N t N)) has exactly the same gaps as a(H per ). 

5. Final remarks 

• The zig-zag nano-tube case has been thoroughly studied by a different 
method in [SB], [27] • The methods employed in this paper seem to be 
significantly simpler than the ones in [2SJ E7J and also apply to all 2D 
carbon nano-structures such as graphene and any single-wall nano-tube. 

After the paper was accepted for publication, the authors received 
the preprint [6j, where the case of armchair nano-tubes is considered by 
methods analogous to the ones of [26], [27] . 

• Other spectral properties of the operators H and H p , e.g. asymptotics of 
gaps lengths or properties of (and formulas for) the density of states can 
be easily derived from the explicit dispersion relations that we obtained 



PETER KUCHMENT AND OLAF POST 



and the well studied properties of the Hill discriminant. As an example, 
we provide a theorem below that describes the smoothness of the poten- 
tial in terms of the gap decay. In order to do so, we call the gaps arising 
as .D _1 ([— |a(p), §a(p)]) the odd gaps G 2 k-i, k — 1, 2, ... and the gaps of 
the Hill operator the even gaps G2k,k = 1,2,.... Notice that we count 
the gaps even when they close (e.g., all odd gaps close for graphene and 
for nano-tubes with integer (pi — p 2 )/3). Let also 7^ be the lengths of 
the gap G k . 

In the theorem below, the operator is either the graphene operator H, 
or the nano-tube operator H p . 

Theorem 5.1. (i) The periodized ID potential q is infinitely differ- 
entiable if and only if ~/ 2 k decays faster than any power of k when 
k — > 00. 

(ii) The periodized ID potential qo is analytic if and only if j 2 k decays 
exponentially with k. 

Since the even gaps are exactly the spectral gaps for the Hill operator 
with the periodized qo, this is an immediate corollary of the results of 
this text and known results of the same nature for the Hill operator 
[T2"l |2"2"| [231 EH EZ] ■ Statements similar to this theorem can be derived as 
easily for other functional classes of potentials, using the corresponding 
results for the Hill operator in [12j. 

Albeit for graphene (as well as for the nano-tubes with p\ — p 2 divisible 
by 3) the absolutely continuous spectrum coincides with the one of the 
periodic Hill operator as a set, the structure of the spectrum is different, 
due to the appearence of the conical singularities inside of each band of 
the Hill spectrum. 

Such singularities can also appear when the even gaps close, but this 
situation is non-generic with respect to the potential q . However, as 
we discussed above, closing the odd gaps and thus appearence of conical 
singularities there is mandatory for any potential in the graphene case, 
as well for nano-tubes T p with p\ — p 2 divisible by 3. 
As we have indicated in the beginning, quantum graphs (quantum net- 
works) have been used to model spectra of molecules at least since 
[5"Tl 153] . However, the validity of such models is still under investiga- 
tion, see e.g., [131 [TBJ [331 SSI HZ1 E2] and references therein. 
The graphene operator H provides also an interesting example in terms 
of Liouville type theorems. As it was established in [38] (see also [37J for 
related considerations), the Liouville theorem for periodic operators holds 
if and only if the Fermi surface consists of finitely many points. Albeit 
this normally occurs at the spectral edges only, it was indicated in [371138] 
that in principle this can happen inside the spectrum. The graphene 
operator provides just such an example. Namely, a direct corollary of 
our results and the ones of |38j is the following 
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Theorem 5.2. Let H be the graphene operator. Suppose D(\) = and 
n > 0. Then the space of solutions u of the equation 

Hu-Xu = 

such that 

Hx)\<c u (i + \x\r 

is finite dimensional. 
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